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\^ , Neutrino propagation in protoneutron stars requires the knowledge 

i of the composition as well as the dynamical response function of dense 

hadronic matter. Matter at very high densities is probably composed of 
other particles than nucleons and little is known on the Fermi liquid prop- 
erties of hadronic multicomponent systems. We will discuss the effects that 
O \ the presence of A hyperons might have on the response and, in particular, 

on its influence on the thermodynamical stability of the system and the 
mean free path of neutrinos in dense matter. 



PACS numbers: 13.15.+g,13.75.Ev,21.30.Fe,25.30.Pt,26.60.+c 
?H \ Core collapse supernovae appear at the last stage of the evolution of 

massive stars PQ. When the nuclear fuel of the star is exhausted, no more 
energy is produced in the interior and the most internal iron shell collapses. 
When the central density reaches the saturation density of nuclear matter 
po, the collapsing matter bounces and forms a protoneutron star (PNS). 
This protostar is expected to reach densities of the order of 2 — 3po anci 
temperatures up to 50 MeV 0. The neutrinos produced during the col- 
lapse are trapped in the dense core. In addition, the entropy increases and 
thus thermal production of neutrino pairs is enhanced. The PNS is then a 
neutrino-rich environment, where neutrinos have been dynamically trapped. 
However, neutrinos can diffuse out to the crust of the PNS, carrying the heat 
of the inner core with them. The early cooling process of the PNS is thus 
driven by neutrinos and it is very sensitive to neutrino interactions in dense 
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matter. For instance, the details of the nucleon-nucleon (NN) interaction 
and the dynamical correlations in the nuclear medium are of a capital im- 
portance for the cooling process [3] . The presence of different particles apart 
from nucleons (such as hyperons) can also influence the neutrino emissivity 
of the star, either because of the effects that hyperons produce on the in- 
medium nucleons or because of the neutrino interactions with hyperons. In 
any case, an important and illustrative measure, which might be useful for 
early cooling process computations or supernovae evolution simulations is 
the neutrino mean free path (MFP) in the hot and dense hadronic medium 
of the PNS. The influence of the different baryons in this quantity has been 
previously studied, mainly within relativistic mean field approaches 

From a theoretical point of view, Landau theory of Fermi liquids ex- 
plains in a successful way the properties of quantum liquids like 3 He or 
nuclear matter [3]. Within Fermi Liquid Theory (FLT), the study of the 
thermodynamical (TD) instabilities or the collective modes of hyper-nuclear 
systems could be easily assessed. The main difficulty faced in the model- 
ing of dense matter lies on our lack of knowledge of the hyperon-nucleon 
(YN) and hyperon-hyperon (YY) interactions as well as on the in-medium 
modifications of baryons. 

The main idea underlying FLT simply states that, under certain assump- 
tions, the properties of Fermi liquids are only related to the excitations of 
the system close to the Fermi surface. One can easily study these excita- 
tions (usually called quasiparticles) and find out the static and dynamic 
properties of the system. In FLT, a fluctuation dn^ = — of the Fermi 
surface will modify the energy Eq of the system in its ground state by an 
amount 5E = E[nk] — E[n^\, which is related with the fluctuation through: 

5E = e{h)6n(h) + i £ fe 2 )M&)M&) , (1) 

One can then define the quasiparticle energy, e(k), as a first variation of the 
energy with respect to the momentum distribution, while a second variation 
leads to identify the function /(fei,/^) to an interaction between quasipar- 
ticles: 

r) 2 F 

f(k u k 2 ) = =, (2) 

5n(ki)Sn[k2) 

Since we are dealing with excitations around the Fermi surface, the quasi- 
particle interaction is usually taken at k\ = k 2 = kp. The only remaining 
dependence is the one in the Landau angle 9±2 between ki and k 2 , which 
can be accounted for by a Legendre expansion: 

/(0i 2 ) = E/^[ cos (M]. (3) 
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The Legendre coefficients Jl are the so-called Landau parameters. A rele- 
vant issue, already pointed out by Landau in his original derivation of FLT, 
is that these parameters can be related to static and dynamical properties 
of the system, such as the effective mass, the compressibility or the spin 
susceptibility. 

Up to now we have restricted ourselves to a system of spinless fermions. 
When other quantum numbers, such as spin are taken into account, the 
quasiparticle interaction should be properly generalized. Usually, one ex- 
pands the quasiparticle interaction in the following operational form: 



Since 1 and a\ ■ c?2 are the S = and S = 1 projector operator in the 
particle-hole channel, it is clear that / and g account for the S = and 
S = 1 spin channels of the particle-hole interaction. Finally, whenever we 
deal with a system of different species of fermions (such as hadronic matter), 
we will have to take into account the presence of the different Fermi surfaces 
and consider a Landau parameter for each pair of species. In the case of 
a system of neutrons, protons and A's, the possible / parameters will be 

jnn jnp jpp f AA 

As a first step in the application of FLT to hadronic systems, we will 
consider a simple case where the in-medium hadronic interaction is of the 
Skyrme type. The NN, NY and YY in-medium interactions will then be 
given by different sets of Skyrme-type parameterizations. Skyrme forces 
are very well known in the NN sector, where they have implied a major 
breakthrough in the study of both the static properties and the collective 
motion of finite nuclei. In the fitting procedure of its parameters, one usu- 
ally includes the properties of nuclear matter and some given properties of 
medium and heavy spherical nuclei. On the other hand, the available exper- 
imental data on heavy hyper-nuclei has triggered the extension of this kind 
of forces to the YN and YY sector and, in particular, several parameteriza- 
tions exist for the AN interaction These Skyrme functionals have been 
usually fitted to the single-particle binding energies of A's in hyper-nuclei. 
To our knowledge, there is only one AA Skyrme-type interaction, fitted to 
the AA binding in double hyper-nuclei [7j. For convenience, from here on 
we will restrict ourselves to a single NN parameterization, the SLylO |Hj set 
of the Lyon group. In the NY sector, we will use the LN1 parameterization 
of Ref. U3, while we will use the different YY parameterizations of Ref. [7j. 

We are interested in the description of hadronic matter in astrophys- 
ical environments. In a neutron star, the environment is defined by the 
/^-equilibrium conditions, which arise from the assumption that the weak 
reactions: 
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Fig. 1. /3-stable composition of a PNS formed of neutrons, protons and A's as a 
function of the total baryonic density. The four panels correspond to two different 
temperatures, T = 20 MeV (left panels) and T = 40 MeV (right panels), and two 
different leptonic fractions, Yi =0.1 (upper panels) and Yi = 0.3 (lower panels). 



p + e «-> A + z/ e (6) 

are in equilibrium in the hadronic medium. In addition, a PNS is a charge 
neutral system (p p = p e ). In Fig. 1 we show some examples of the /3-stable 
composition of protoneutron stars at T = 20 MeV (left panels) and T = 40 
MeV (right panels) temperatures as a function of density. The upper (lower) 
panels correspond to a leptonic fraction of Y[ = 0.1 (Yi = 0.3). Different sets 
of NN and YN forces, however, would produce different particle fractions 
and different astrophysical scenarios. 

As an example of the application of FLT to hadronic systems, we will 
compute the response function of hadronic matter from the Landau param- 
eters. If we consider diffusion reactions, in which neutrinos create a neutral 
particle-hole excitation on the system, the lowest-order response function is 
given in linear response theory by a 3 x 3 matrix structure: 

/ lift" \ 
n = Eg" , (7) 

V o o n£ A / 

where each diagonal ITq is given by the corresponding Lindhard function: 
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Thermodynamical stability 
T=20 MeV, Y,=0.3, SLylO + LN1 




Fig. 2. Stability criteria for hadronic matter for the spin S — (right panel) 
and S — 1 channels (left panel) as a function of the total baryonic density. As 
illustrative examples, we have chosen a temperature and leptonic fraction of T = 20 
MeV and Yi = 0.3, respectively. Different sets of AA interactions are shown with 
solid (LLO), dotted (LL1), dashed (LL2) and dash-dotted lines (LL3). 



The long range correlations are included through the Bethe-Salpeter equa- 
tion: 



n = n + n K es n . 



(9) 



This is the well-known RPA approximation. The residual particle-hole inter- 
action V res is obtained from Eq. © and reduced to the monopolar Landau 
parameters / (g) in the S = (S = 1) channels. The Bethe-Salpeter equa- 
tion becomes then a coupled set of algebraic equations where the residual 
interaction is given by a 3 x 3 symmetric matrix. In the S = channel, one 
finds: 

/ inn 
I JO 

f pn 

JO 

V fo An 

The analytical expressions for the different / and g parameters can be found 
from the second derivatives of the energy Eq. @ and are explicitely given 
0. The solution of the Bethe-Salpeter equation is then: 




m 



n = (i - u v r( 



(10) 



If the determinant of the matrix (1 — IloV^ es ) becomes singular, the equa- 
tion might not be solvable. In the low energy limit, one can see that this 
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determinant reduces to the determinant of the TD stability matrix: 



det [1 - V res U 
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A change of sign in this determinant signals the onset of a TD instability. 
Thus, the Landau parameters of hadronic matter are related to its TD 
stability. 

Results for the values of the low-energy limits of these determinants are 
shown in Fig. 2, where they are plotted as a function of the total bary- 
onic density in /3-stable conditions. The different sets of AA interactions 
give different determinants due to both the changes in the /3-stable particle 
fractions and the differences on the parameters of the force. As a result, dif- 
ferent TD instabilities arise depending on the YY interaction. In particular, 
the onset of a ferromagnetic transition (signaled by the zero of det[l + Go]) 
strongly depends on the set of forces that is used. In this case, all the 
forces tend to a zero value at high densities (all of them are thus close to 
a ferromagnetic instability), but only LL3 suffers from a real ferromagnetic 
transition at p ~ 0.74 fm~ 3 . 

Once the Bethe-Salpeter equation is solved, we can find the structure 
functions by means of the finite temperature relation (obtained from the 
detailed balance relationship): 



1 -exp[-w/T] 
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where (S) denotes the two possible spin channels, S = and S = 1. In the 
calculation of the MFP we shall use the following "generalized" structure 
function, where all the channels are taken into account by the following 
matrix product: 
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(12) 



and where the are the corresponding weak charges of the particles which 
can be found in 

In the left (right) panels of Fig. 3 we show the structure function for the 
S = (S = 1) channel as a function of the energy, for a fixed momentum 
of q = 30 MeV. For the spin dependent channel, a collective mode around 
u ~ is observed when the density increases: this is again a signal of a 
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RPA structure functions 
T=20 MeV, Y =0.3, SLy 10 + LN 1 + LL2 



— p=0.16fm" 3 




l ji 


p=0.32 fm" 3 




I ii 


" — p=0.64fm" 3 




'! i' 


- - p=0.96fm" 3 














-20 
CO [MeV] 



-20 20 
CO [MeV] 



40 



- 300 B 



3 

- 100 



Fig. 3. Structure functions for the S=0 (right panel) and S=l (left panel) channels 
as a function of energy for a transferred momentum of q — 30 MeV. The calculations 
have been performed in the same /3-stable conditions as for the previous figure at 
four different densities. 



close-by ferromagnetic transition. One can see that when such a strong col- 
lective mode dominates the structure function, the neutrino MFP is strongly 
suppressed. Thus, a different choice of NL and LL interactions which mod- 
ifies the TD instabilities of the system, might have a great influence on the 
neutrino MFP. 

Finally, the mean free path is given by the following integral jlUj : 



A ~ 4^2 



/ dw-f dqq[l-n(El)} 

J oo J-Jv J I Id I 

1 -(- cos 0)<S( O ) (a>, (?) + (3 — cos0)<S(!)(u;, q) . (13) 



We show in Fig. 4 the results for the mean free path of neutrinos in dense 
matter for a leptonic fraction of Yi = 0.3 and two temperatures, T = 20 
MeV and T = 40 MeV. The left panels correspond to the mean-field MFP 
computed from the non-correlated response functions ITg . The solid lines 
correspond to /3-stable matter composed by nucleons, electrons and A's, 
while the dashed lines correspond to matter without hyperons. The neu- 
trino energy is E v = 3T MeV. For both temperatures, a clear reduction 
of the MFP is observed when A's are included in the system. This is due 
to the opening of a new channel of interaction, with which neutrinos can 
interact, thus increasing their cross section (and reducing the corresponding 
MFP). The right panel of Fig. 4, on the other hand, illustrates the ratio be- 
tween the MFP computed with RPA correlations and the mean field MFP. 
The density ranges explored for each MFP correspond to the ranges where 
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Fig. 4. Left panel: neutrino MFP in the mean field approximation for /3-stable 
matter with Y\ = 0.3. The full (dashed) lines correspond to matter with (without) 
A's. Right panel: ratio of the RPA and mean field approximations of the neutrino 
MFP for the same conditions. 



no instabilities are observed for the corresponding /3-stable matter. A spin 
instability is observed for the matter formed by nucleons and electrons at 
densities above p ~ 0.5 fm~ 3 . The presence of A's pushes the ferromag- 
netic instability to higher densities, and thus a non-zero neutrino MFP is 
observed in the high density range. This kind of effects, where the presence 
of A's in the system result in a suppression of the instabilities of matter, 
are probably the most relevant ones for the neutrino propagation in hot and 
dense hadronic matter. 
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